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Abstract 



CLl 



We explore the possibilities for measuring the quark content of the proton and neu- 
tron using neutrino beams produced at a muon storage ring. Because of the nature 
^h. of the beams, small nuclear targets such as hydrogen and deuterium can be con- 

sidered, as well as polarized targets. The statistics expected from these targets are 
calculated using nominal muon storage ring intensities, and the resulting statistical 
errors on the numerous structure functions available are given, for both polarized 
and unpolarized targets. It is shown that with a relatively small target, the struc- 
ture functions F2, X-F3, xg\ and xg^ for neutrinos and antineutrinos on protons 
and deuterium, either unpolarized or polarized, could be determined with excellent 
precision over most of the accessible kinematic range. 



X 



1 Introduction 



Deep Inelastic Scattering (DIS) has long been the definitive process for the 
determination of the quark content of protons and neutrons. Charged lepton 
scattering has dominated the field in terms of precision determinations of 
the sum of the quark and antiquark distributions, and the associated gluon 
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Fig. 1. The kinematic region in the (x,Q 2 ) plane available at a 50 GeV neutrino 
factory. 

distribution, but neutrino scattering has thus far contributed complementary 
measurements of the valence quark distributions, as well as measurements of 
the strange sea. However, neutrino scattering has in the past been plagued by 
tiny interaction rates and large beam spot sizes, requiring targets on the order 
of several meters wide and several hundred tons to get appreciable statistics. 



With the advent of a muon storage ring the flux of neutrinos at a near detec- 
tor would be several orders of magnitude higher than at present experiments, 
and concentrated in a much smaller spot size. Because of this one can now 
consider using compact hydrogen and deuterium targets, rather than iron. 
These targets have the advantage of allowing measurements of the valence 
quark distributions without nuclear effects, or conversely one can finally mea- 
sure nuclear effects in valence quark distributions by comparing results using 
different targets. Many of these ideas (as well as other high-rate neutrino ex- 
periments at muon storage rings) are considered in references [3], [4], and [5]. 
Because it is expected that the storage ring will run in roughly equal running 
times in fi + and /i~ mode, the fluxes for v e and v e will be approximately equal, 
as will the fluxes for v^ and v^. In conventional beam neutrino experiments 
the dominant statistical error has been the antineutrino event rate, because 
the typical total antineutrino event rate (on the targets used) has been only 
20-25% of the neutrino event rate. 



Aside from just statistical considerations, the neutrino beams from a muon 
storage ring offer ways to lower previously important systematic uncertainties 
as well. Because of the well-known incoming v spectrum one has an accu- 
rate determination of the beam energy, as well as a calibration tool for the 
detector. One also has an extremely pure neutrino beam in terms of sign se- 
lection; if /i~(/x + ) are circulating in the ring the fluxes available are v^ and v e 
(z/ M and v e ). By identifying the flavor of the final state lepton in a charged 
current interaction one knows if the initial lepton was a neutrino or an an- 
tineutrino. Finally, again because of the small beam spot size, polarized Z/-DIS 
experiments could be performed for the first time, allowing a full flavor de- 
composition of the nucleon spin. Because these measurements would still be 
inclusive DIS measurements, they would not suffer from the fragmentation 
uncertainties intrinsic to semi-inclusive measurements such as those made at 
HERMES. 

One disadvantage of the neutrino beams from a muon storage ring is that 
because the rings are expected to operate at a relatively low beam energy 
(30 to 50 GeV), a lower and smaller range of momentum transfers will be 
available. The expected kinematic range for a 50 GeV beam is shown in figl. 

The remainder of this paper consists of a description of the theoretical frame- 
work by which the quark structure of the nucleon can be measured, a brief 
discussion of what the neutrino target and detector would look like, and finally, 
preliminary leading order estimates of expected uncertainties for a variety of 
stucture functions. 



2 Unpolarized Structure Functions 



Unpolarized charged current structure functions are defined through the de- 
composition of unpolarized differential charged current cross-sections into in- 
variant functions of the momentum of the struck quark (x) and and the mo- 
mentum transfer squared of the W boson (Q 2 ): the standard definitions give 



where S = 2mE is the centre-of-mass energy, E is the neutrino beam energy, 
assumed to be ^> m, and the ± signs refer to the sign of the charged cur- 
rent: W + exchange for v scattering and W~ for V. y is the fractional lepton 
energy loss, or {E v — Eg)/E u . In neutrino scattering, x,y, and Q 2 can all be 
determined simply by measuring the outgoing lepton energy and direction, 
and the hadronic energy in the event. There are then in principle six indepen- 
dent structure functions to be measured for every target. The proton structure 



functions are the same as those measured at higher Q 2 (see fig.l in charged 
current e ± p scattering at HERA. 

In the parton model four of the these structure functions are related through 
the Callan-Gross relations F 2 = 2xF\. the longitudinal structure function 
Fl = F 2 — 2xFi begins at 0(a s ) in perturbation theory. The six structure 
functions Fi' u , F% ,v and F$' v may in the parton model be expressed in terms 
of parton densities as 



F? = u + d + s + c, F? = u + d + s + c, 

F 2 u = 2x{u + d + s + c), Fl = 2x{u + d + s + c), (2) 

xF% = 2x{—u + d + s — c), xF£ = 2x{u — d — s + c), 

where we have set the CKM mixing angles to zero for simplicity, and restricted 
attention to the first four flavors. To go from a proton to a neutron target 
(assuming isospin invariance) we interchange u and d. It is thus not difficult 
to see that by constructing appropriate linear combinations of all eight inde- 
pendent structure functions (conventionally taken as {F2 V ) p , n and {xF^' v ) p>n ) 
obtained by v and v scattering on proton and neutron (or deuteron) targets it 
is possible to separately disentangle u±u, d±d and s±s provided only that we 
can determine c±c), either theoretically, or else empirically by tagging charm 
in neutral current processes (as is done currently at HERA). More explicitly, 
assuming isospin invariance between proton and neutron targets we have 

(FZ +P ) P = {FZ+% = x(u + u + d + d + s + s + c + c), 
(xF^-% - (xF^-% = -2x(u + u-(d + d)) } 
(xF»-% + (xF»-% = 2x(s + s-(c + c)), 
(xF^ +p ) p = (xF^ +p ) n = x{u - u + d - d + s - s + c - c), (3) 

(F»~% - (F»-% = -2x(u -u-(d- rf)), 

(Fr% + (F»-% = 2x(s -s-(c- c))), 

where F^ ±u = \{F" ± F"). The first and fourth of these equations are the 
structure functions F£ +u and F% +v normally measured in neutrino scattering 
(though on heavy targets), while the second and third allow flavour decom- 
position of the total q + q distributions, and the fifth and sixth a similar 
decomposition for the valence distributions. To separate out strangeness from 
intrinsic charm empirically would require either a tagging of charm in the fi- 
nal state to give an independent determination of strangeness alone, or (in 
principle at least) a combined analysis with neutral current structure function 
data. 



2 It is inappropriate to consider only three flavors in charge current scattering since 
scattering off a strange quark produces a charmed quark in the final state. 
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Fig. 2. A strange quark distribution extracted from charged current cross section 
data [6]: results from a NLO CCFR dimuon determination [7] are shown for com- 
parison. 

In perturbative QCD it is possible to perform a complete NLO global analysis 
using well known results for NLO coefficient functions and anomalous dimen- 
sions. The charm contribution can then be computed perturbatively (on the 
assumption that the intrinsic charm is very small). Since only F% +v contains 
a singlet component, while F^~ v and F^ ±u are entirely nonsinglet, a clean 
extraction of both a s [1] and the gluon distribution from scaling violations is 
possible with data taken over a sufficiently wide kinematic range. To perform 
such an analysis it is necessary to either extract the structure functions from 
the differential cross section in a model independent way [2], or simply to fit 
directly to the cross section [6]. Fits to BEBC,CDHS and CDHSW data sug- 
gest that it is already possible to extract s + s (see fig. 2) and possibly even 
s — s: it will be interesting to see the results of this type of analysis applied 
to the much more precise CCFR/NuTeV data. 



3 Polarized Structure Functions 



Polarized structure functions may be defined in analogy with the unpolarized 
ones through asymmetries in the polarized cross- sect ions: for longitudinal po- 
larization we may write 
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where the polarization asymmetry Aa = aZ, — cx^Z (the unpolarized cross- 
section (1) being a = |(a^ + cr^z), which accounts for the extra factor of 2 on 
the right hand side of (4)). With these definitions [8] 0in the high energy limit 
E ^> m g 2 and g% drop out, and we are left with an expression of the same 
form as the unpolarized decomposition (1), but with iq — > g§, F 2 —>■ g^ and 
F3 — > 1g\. Thus we again have six partonic structure functions for every target 
and for high energy polarization asymmetries the complicated decomposition 
(4) then becomes simply 
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1,1 A " - >±—\±y(i-mx g r + (i-y)gr + y 2 xgn (5) 



dxdy 7r(l+Q 2 /M^ ) 



The remaining four structure functions g 2 v and g^' v have no simple partonic in- 
terpretation and are contaminated by twist three contributions: their twist two 
components are fixed by the Wandzura-Wilczek relation (giving g 2 in terms 
of <7i) and a similar relation [8] which gives #3 in terms of g±. They are most 
easily determined by measuring transverse asymmetries: such measurements 
are very difficult however because the transverse asymmetry is suppressed by 
m/Q. 

In the parton model g^ and g$ are related by an analogue [9] of the Callan- 
Gross relation: g± = 2xg$(l+0{a s )). The flavor decomposition of the structure 
functions g v { v , g% v and g^' u may thus be expressed in terms of parton densities 

as 



g\ = Au + Ad + As + Ac, g% = Au + Ad+ As + Ac, 

g% = 2x(-Au + Ad + As- Ac), g\ = 2x(Au - Ad- As + Ac), (6) 

gl = -Au + Ad + As - Ac, g\ = Au - A d - As + Ac, 

in precise analogy with the unpolarized case (3): comparing (1) with (5), 
F\ ~~ * 95, F2 -^ 94, |-^3 - > 9\ and q — > Aq, q — > —Aq (changing a quark to 
an antiquark also flips its helicity). Again, by constructing appropriate lin- 
ear combinations of all eight independent structure functions (conventionally 
taken as (gi ,l %, n and (g^' u ) Pjn ) obtained by longitudinally polarized v and v 
scattering on proton and neutron (or deuteron) targets it is possible to sepa- 
rately disentangle Am ± Au, Ad ± Ad and (As ± As) just as in eq.(4). 

Some combinations of the polarized structure functions are of particular in- 
terest. For example, writing g\ = \{g" ± gX), the first moment of 

2{g v l + % = 2{gl +9 ) n = Au + Au + Ad + Ad+ As + As + Ac + Ac (7) 



3 There are many variants in the literature: see [8] for a compilation. 



is the axial singlet charge ao- This is a much more direct measurement than 
the traditional one through electron-proton or deuteron DIS since in the latter 
case one must first subtract the octet charge a§ which is then only determined 
indirectly through hyperon decays (see [10] for a recent review). Thus in Z/-DIS 
one would have a direct check on the anomalous suppression of ao- Similarly 
first moments of 



{9V)v-^V)n- 


= -(Au + Au) + Ad + Ad, 


(gr)p + (gr%-- 


= As + As- (Ac + Ac), 



give direct measurements of the axial charge 0,3 (again currently only measured 
indirectly through /5-decay) and of the contribution of strange quarks to the 
nucleon spin, as would the tagging of charm in the final state. Flipping the 
signs, we can also determine the contribution of valence quarks to the spin, 
since 

2(gl +9 ) p = 2(^ +p ) n = Am - Au + Ad - Ad + As - As + Ac - Ac, 
{gl-% - (g»-% = - ( Au - Au) + Ad - Ad, (9) 

{9i~% + (9i~% = As - As- - (Ac - Ac), 

so one could even check for intrinsic strange polarization As — As. None of 
these valence polarizations can be cleanly measured in current polarization 
experiments. 

In practice this flavor separation would be best performed by a global fit in 
NLO perturbative QCD: all the NLO anomalous dimensions [11] and coeffi- 
cient functions [12-14] are known, the latter for heavy quarks, so the polarized 
charm contribution can be computed perturbatively. 



4 Estimates of Neutrino interactions and Event Rates 



In order to estimate the precision with which the structure functions above 
can be measured at a neutrino factory on nuclear targets, a GEANT-based 
Monte Carlo was used, with hydrogen and deuterium chosen as targets. This 
simulation used LEPTO and JETSET (versions 6.5.1 and 7.408 respectively) 
to determine what the struck quark and its fragmentation were for each inter- 
action, and the parton densities used were from CTEQ2MS [15]. The hadron 
showers were not traced through the detector; because of the low densities of 
the targets and their sizes it is expected that most of the shower will leave 
the target and be collected in a low mass particle-tracking system. Therefore, 
in the results described below there is no smearing from detector resolution. 
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Fig. 3. Number of muon neutrino events per GeV per g/cm 2 for 10 20 muon decays 
in an 800m long straight section followed by a 40m shielding section and a detector 
of radius 5 or 50 cm. Other storage ring parameters are given in the text. 

The acceptance for muon neutrino and antineutrino charged current inter- 
actions is assumed to be 100% for events with muons above 3 GeV, and to 
remove events from the quasielastic and resonance region only events with 
hadron energies above 1 GeV were considered. It should be noted that the 
precisions listed here are pessimistic, since they are calculated assuming only 
an incoming ^(P^) flux. The P e (z/ e ) fluxes that arrive simultaneously could 
in principle donate roughly another factor of two in event statistics. However, 
since the acceptance and backgrounds for v e charged current events are much 
more detector dependent these events are not considered here. 

The question of where the high-rate neutrino physics experiments would occur 
at a neutrino factory is not a trivial one, but for a rough estimate of the overall 
statistics in neutrino (fi~ in the storage ring) and antineutrino (/i + in the 
storage ring) running, consider two scenarios: a 25 and 50 GeV muon storage 
ring with 800m straight sections, followed by 30m of active shielding. Two 
target sizes were considered, namely a 5cm radius or a 50cm radius. In both 
cases the muon beam spot size was 1.23cm(x) x 0.883cm(y), and the muon 
beam divergence was 0.73 mrad in both the vertical and horizontal directions. 
Figure 3 shows the muon neutrino fluxes for these different scenarios. The 
(x,Q 2 ) regions accessible are comparable between two different radii at the 
same energy: the most important factor is the loss of statistics for the smaller 
target. Note that the rate difference between the two radii is roughly 10 at 
high energies rather than a factor of 100. Table 1 gives the neutrino charged 



Muon Energy 


Detector Radius 
5cm 50cm 


25 GeV 
50 GeV 


24.5K 841K 
131K 2900K 



Table 1 

Muon neutrino charged current interaction rates for 10 20 muon decays (one year) 

for detectors of different radius and for storage rings of different energies. 

current interaction rates for these different scenarios on a 1 g/cm 2 isoscalar 
target. 

Given that the density of cryogenic deuterium is at least 0.162 g/cm 3 , 1 g/cm 2 
corresponds to only 6 cm of deuterium. Targets 1.3 meters long of cryogenic 
(polarized!) material have already been used by SMC, so one can imagine 
multiplying the statistics listed in table 1 by a factor of 20 and still have 
a 10g/cm 2 target after fiducial cuts on the vertex. Much larger targets of 
liquid Hydrogen have been used in the past in neutrino experiments, but with 
extremely low neutrino event statistics. 

The relative cross sections per nucleon between vD 2 : vD<i '. vH 2 '■ vH 2 are 
approximately 2:1:1.3:1.3, (as calculated by GEANT) due to the higher abun- 
dance of up quarks in H 2 compared to D 2 . For an unpolarized target one 
can consider radii of 50cm. Due to the strong B field and low temperature 
requirements for the polarized targets (as well as the small beam spot size), 
at present only polarized targets of 5cm diameter have been used for charged 
lepton scattering experiments [16]. Given the rate of advances in cryogenic and 
magnetic field technology, it is not unreasonable to expect that much larger 
targets will be available several years from now. For that reason we consider 
targets of polarized materials that are the same size as the modest unpolarized 
targets. 



5 Unpolarized Cross Sections at a Neutrino Factory 



To estimate the uncertainties on the unpolarized structure functions at a muon 
storage ring experiment, a leading order "model-independent" analysis was 
performed on the Monte Carlo data, as is described in [2]. Simply put, the 
differential cross sections are measured in (x, y, E) bins, and then they are 
collected into (x, y, Q 2 ) bins, such that for a given (x, Q 2 ) bin there is a range 
of neutrino energies contributing as y goes from to 1. The errors on the cross 
sections in each bin are taken to be the statistical errors alone resulting from 
the above fluxes (taking into account the total cross section differences between 
the different targets and probes). For a given (x, Q 2 ) bin the y distribution is 
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Fig. 4. Expected results for F 2 ' u from one year of running in each mode, with a 0.1 
g/cm 2 target: for a 10 g/cm 2 target the errors would be reduced by a factor of 10. 

fit to the function in eqn.(5), setting 2xFi = F2. This will modify the central 
values fit since in reality the ratio R = (F% — 2xF\)/ F% is non-zero in neutrino 
scattering just as it is in electron scattering [17], but the errors on F 2 and xF 3 
should not be affected. Also, there are no bin-centering corrections applied 
here, but the resulting fit errors should again not be affected. 

Figures 4 and 5 show the fit results for the structure functions F 2 and xF 3 on 
both deuterium and hydrogen for both neutrinos and antineutrinos. Because 
of limited Monte Carlo statistics, the errors shown above are for 0.1 g/cm 2 



10 




Fig. 5. The same as fig. 4, but this time showing results for xF%. 

targets, for 10 20 50 GeV muons decaying in the straight section in each mode 
(/i + and /j,~), and for a target of radius 50 cm. The \ 2 for the fits were roughly 
at 1 per degree of freedom or better, because the Monte Carlo statistics were 
slightly higher than the expected experimental statistics for these targets. 
The error reported by the fit does represent the experimental error expected, 
however, not the error due to the Monte Carlo statistics. 

These data access a Q 2 range that is lower than most of the CCFR data. For 
the expected precision of a modest-sized target, one would divide these errors 
by a factor of 10 which assumes a 10 g/cm 2 fiducial target. What makes these 
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Fig. 6. Estimated errors on F^ v and xF% u after one year of running on each 
of neutrinos and antineutrinos, on a target of 10g/cm 2 of D 2 and H 2 — D 2 for the 
range of x and Q 2 space above Q 2 = 1 GeV . Note that one can achieve errors of 
better than 0.01 on sums and differences of structure functions with these statistics 
in most of the kinematic range accessed. 



measurements unique however is the possibility of taking particular sums and 
differences between combinations of v and v, hydrogen and deuterium struc- 
ture functions. Figure 6 shows the expected errors on the various structure 
function combinations listed in eqn.4, remembering that D 2 = |(p + n), while 
(H 2 -D 2 ) = ±(p-n). 
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The errors on the sums and differences of F 2 and xF% are comparable, but 
while the sums are of order unity, some of the differences are expected to be 
quite small. Also, the errors on xF 3 are larger than those for F 2 for every 
combination, although they approach each other as Q 2 increases. What is 
important here however is that the errors are uniformly small; less than 0.01 
for most of the (x, Q 2 ) range accessible. For comparison, the statistical errors 
on (u + i?) Fe structure functions from CCFR in its kinematic range vary from 
0.01 to 0.04 [1]. 



6 Polarized Cross Sections at a Neutrino Factory 



The polarized cross section measurements will be more difficult due to sev- 
eral factors, Firstly, to keep the target polarized it must be in both a strong 
magnetic field and a very low temperature container, so the target size may 
be limited to smaller volumes than for unpolarized targets. Secondly, since 
one is taking differences of cross sections the polarized cross-sections will be 
smaller than the unpolarized ones and the fractional error correspondingly 
larger. Finally, the fraction of polarized nuclei will necessarily be lower than 
100%. Two possible polarized target materials used in the past in charged 
lepton scattering are polarized solid butanol (used by the SMC collaboration) 
[16] or a polarized HD target (used by the LEGS collaboration) [18]. 

If one only polarizes the hydrogen or the deuterium in the target sample 
then when one takes the differences between opposite polarizations the result 
will only depend on the polarized nuclei. In this way one can use deuterated 
butanol and butanol to measure the D and the H cross section differences. 
In the HD target, the two components are independently polarizable. So for 
measurements of p+n or D alone, one can simply polarize the D in the sample 
and leave the H unpolarized. For measurements of p-n, one can either polarize 
only the H, or use targets where the the H and the D are polarized in opposite 
directions. In the latter case one is effectively scattering off polarized neutrons. 

To see how the errors on the polarized cross sections compare to those on the 
unpolarized cross sections, consider the following argument: for a perfectly 
polarized H 2 target of the same size as that considered above (10g/cm 2 , 50cm 
radius) if all the protons were polarized, then the error on the cross section 
difference would rougly be l/\/2 times the error on the unpolarized case, as- 
suming one integrated 10 20 muon decays in one polarization and 10 20 muon 
decays in the opposite polarization. Factors such as the incomplete polariza- 
tion of the target and the fact that the polarized target is not 100% H 2 enter 
into the error in the following way: 
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Characteristic 


p-Butanol 


D-butanol 


HD(D^) 


HD (HftD^) 


Density (g/cm 3 ) 


0.61 


0.69 


0.05 


0.05 


D/H Polarization 


86% 


51% 


70% 


70%(D), 
95% (H) 


Dilution Factor v 


0.15 


0.22 


0.4 


1.0 


Dilution Factor v 


0.25 


0.22 


0.55 


1.0 


Length (cm) 


120 


120 


10 


10 


(g/cm 2 ) 


73.2 


82.8 


0.5 


0.5 


Diameter (cm) 


5 


5 


3 


3 


B Field 


2.5 T 


2.5T 


7T 


7T 


Temperature 


0.1K 


0.1K 


1.5K 


1.5K 


total target factor (y) 


2.6 


4.4 


4.3 


2.6 


total target factor [p) 


1.6 


4.4 


3.1 


4.8 



Table 2 

Comparison of four polarized targets: p-butanol, D-butanol, an HD target with only 
the D polarized, or HD with both H and D polarized oppositely. The 'density' is the 
effective density for solid butanol (which depends on the packing fraction), and the 
dilution factors were calculated based on neutrino and antineturino cross sections 
on protons and isoscalar nuclei. The total target factor only takes into account the 
ratio between the given material and liquid H 2 or D 2 densities, the polarization, 
and the dilution factor. 



Cpol 



^"unpol /Punpol 
fu,&Py Ppol 



(10) 



where p U npoi/p P oi is the ratio of target densities to H 2 or D 2 , j v ^ is the dilution 
factor and P is the polarization of the H 2 or D 2 . The dilution factor f VjP is the 
v (or v) cross-section weighted ratio of the polarized nucleon to total nucleon 
content of the target. Because the different targets have of necessity different 
ratios of protons and neutrons, the neutrino and antineutrino dilution factors 
will be different. 

Table 2 compares the SMC and LEGS targets for both the D 2 cross-section 
difference and the H 2 cross-section difference [16] [18] [19]. In summary, in the 
SMC target the density is high but the dilution factor is small, while for the 
HD target the density is low but the dilution factor is high. The HD target is 
newer and has not been made in as large samples as the SMC target, but the 
HD targets are modular and one can easily imagine adding several together. 
Overall the multiplicative factors coming from target details alone range from 
1.6 to 4.8. 



14 



o 



o 



10 



10 



Q 2 = 1 .25 GeV 2 



x 
o 



10 r 

-2 : 

10 r 



10 r 



o 
c 

25 



o 
o 

c 
o 

i_ 

o 



10 



-2- 



1- 



Q 2 = 2 GeV 2 




Q 2 = 3.16 GeV 2 



10 

-2: 

10 r 

-3: 

10 - 

10 r 

-2: 

10 r 

-3: 
10 r- 



Q 2 = 5 GeV 2 



Q 2 = 7.9 GeV 2 



0.10.20.50.40.50.60.70.80.9 1 

x 

— err xg 5 (p±u) D 2 — err xq^{v±v) D 2 

— err xg 5 (i/± z>) (H 2 — D 2 ) — err xq^(v±V) (H 2 — D 2 ) 

Fig. 7. Error on polarized structure functions on xg\ v and xg§ v coming from 
one year of running for each combination of neutrinos, antineutrinos, polarized solid 
butanol and deuterated butanol in 2 polarizations for one small range of x and Q 2 
space. Note that one can achieve 0.05 errors or better on spin structure functions 
with these statistics. 

To translate between the cross section uncertainty and the polarized structure 
function uncertainty, recall that the functional form (up to terms in m/E 
which are going to be small even for a 50 GeV muon storage ring) for the 
unpolarized cross sections (eqn. 1) is the same as that for the polarized cross 
sections (eqn. 4), up to the substitutions F 2 — > g± = 2xg 5 and F 3 — > 2gi. 
Therefore, the error on xgi or xg 5 for a given x, Q 2 bin is the same as that of 
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xF^/2 or F2/2 in the unpolarized case multiplied by the ratio of the errors in 
the cross sections and the extra \/2 from having two targets, described above. 

Figure 7 shows the resulting errors on the various linear combinations of polar- 
ized structure functions, listed in eqn.4, remembering again that D 2 = |(p+n), 
while H2 — D2 = \{p — n). We have assumed one year's running at each po- 
larization, each neutrino flavor, and each target. Note that at the very least 
it is assumed that one would run in two target polarizations simultaneously, 
so this represents a total of 4 to 8 years of running at 10 20 /i ± decays per year. 
The targets are the same size as in figure 6 but filled with either butanol or 
deuterated butanol, as was done by SMC [16]. 

The errors on the sums and differences of xg\ and xg§ are comparable, but 
while the sums are of order unity, some of the differences are again expected 
to be quite small. Also, the errors on xg§ are smaller than those for xg\ for 
every combination, although they approach each other as Q 2 increases. What 
is important is that the errors are uniformly small, just as in the unpolarized 
case: less than a few times 0.01 for most of the accessible (x, Q 2 ) range. 



7 Conclusions 



We have examined the novel structure function measurements that could be 
made at a neutrino factory, and provided estimates of their precision as a 
function of x and Q 2 for nominal storage ring running (10 20 muons per year in 
a 50 GeV muon storage ring). Because of the clean separation between valence 
and sea afforded by neutrino and antineutrino running and the possibility of 
using both deuterium and hydrogen targets, such experiments could at last 
determine flavor by flavor the valence and sea quark distribution functions with 
statistical errors uniformly of order 0.01 in each bin. Systematic errors due to 
nuclear effects and beam energy would be minimal. Furthermore, by running 
with polarized targets already developed by the charged lepton scattering 
community the spin components of the proton, quark by quark, could be 
determined for the first time, again with statistical errors uniformly of order 
a few times 0.01. Such measurements would resolve definitively the questions 
raised by the EMC experiment, in particular by determining both the strange 
and gluon contributions simultaneously in one experiment. Combining this 
data with complementary data from eRHIC or polarized HERA could result 
in precise tests of our understanding of the spin structure of the nucleon. 

A 50 GeV neutrino factory would be an ideal partonometer, revealing the 
partonic structure of the nucleon in exquisite detail. We hope that it may be 
possible to enjoy such a facility in the not too distant future. 
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